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By Stuart C. Brown 

SUMMARY 


A method is presented for calculation of static aeroelastic effects 
on wings with supersonic leading edges and streamwise tips. Both chord- 
wise and spanwise deflections are taken into account. Aerodynamic and 
structural forces are introduced in influence coefficient form; the former 
'are developed from linearized supersonic wing theory and the latter are 
assumed to he known from load-deflection tests or theory. 

The predicted effects of flexibility on lateral -control effectiveness, 
damping in roll, and lift-curve slope are shown for a low-aspect-ratio 
wing at Mach numbers of 1.25 and 2.60. The control effectiveness is shown 
for a trailing-edge aileron, a tip aileron, and a slot-deflector spoiler 
located along the 0.70 chord line. The calculations indicate that the 
tip aileron is particularly attractive from an aeroelastic standpoint, 
because the changes in effectiveness with dynamic pressure are small 
compared to the changes in effectiveness of the trailing-edge aileron 
and slot-deflector spoiler. 

The effects of making several simplifying assumptions in the example 
calculations are shown. The use of a modified strip theory to determine 
the aerodynamic influence coefficients gave adequate results only for the 
high Mach number case. Elimination of chordwise bending in the structural 
influence coefficients exaggerated the aeroelastic effects on rolling- 
moment and lift coefficients for both Mach numbers. 


INTRODUCTION 


The characteristic low-aspect-ratio thin wings used for supersonic 
flight have aerodynamic and structural properties which differ consider- 
ably from the high-aspect-ratio subsonic wings for which methods for 
obtaining static aeroelastic effects, such as reference 1, have proved 
to be adequate. When thin wings are subjected to high dynamic pressure, 
there are large aeroelastic effects, which may include chordwise as well 
as spanwise deflections. Consequently, practical methods are needed to 
obtain supersonic aerodynamic influence coefficients and to calculate the 
resulting aeroelastic effects. 
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Several investigations of effects of chordwise deformations on 
aeroelastic effects at supersonic speeds have been made. One approach 
is to express the aerodynamic and structure!. characteristics in equation 
form, combine them, and then solve the combined equation for equilibrium 
conditions (e.g., refs. 2 and 3)- However the equations are only tract- 
able for cases in which the structural and aerodynamic characteristics 
can be expressed in simplified form. An a_.ternate approach, which is 
applicable to more general cases, is to obuain separate sets of structural 
and aerodyn ami c properties in the form of influence coefficients which 
define, respectively, the distribution of deflections and loads due to a 
load or deflection at specified points on a grid network. The resulting 
equations for static equilibrium are then solved numerically. Numerical 
methods have been developed for obtaining unsteady aerodynamic influence 
coefficients for use with flutter analysis (e.g., refs. 4 and 5)- How- 
ever, these methods which were derived for the more mathematically diffi- 
cult problem of determining unsteady aerodynamic forces may be somewhat 
cumbersome when applied to the static case. Consequently, the aerodynamic 
influence coefficients which appear in this report should be useful for 
pseudostatic aeroelastic calculations which occur at frequencies where 
unsteady lift effects are small. 

In the present report, a method, based on linearized potential flow 
theory, for the calculation of static aerodynamic influence coefficients 
is presented for wings with swept supersonic leading edges and streamwise 
tips. Both chordwise and spanwise distributions of loads and deflections 
are included. A numerical iteration method for combining the aerodynamic 
influence coefficients with a set of structural influence coefficients to 
calculate the resulting load distribution over a flexible wing is also 
given. The application of this method is Illustrated by the calculation 
of the rolling performance of several type 5 of lateral-control devices 
for a low-aspect -ratio thin wing with a supersonic leading edge at Mach 
numbers of 1.25 and 2 . 60 . 


SYMBOLS 


[A] aerodynamic influence coefficient natrix 

(An element in the ith row and jth column indicates a 
load at station i due to an angle of attack at station j.) 

[B] matrix for converting to dimensionLess load Coefficients over a 

chordwise interval (eq. (9)) 

lift 

0^ lift coefficient , 


C ^a 


rate of change of lift coefficient with angle of attack, per 
radian 
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„ n . . . _ _ . . , rolling moment 

rolling-moment coefficient, 

C i rate of change of rolling-moment coefficient with wing- tip helix 

P , £b ' 

angle , 

Cp pressure coefficient 

[D] differentiating matrix used in converting to aerodynamic influence 

coefficients (eq. (ll)) 

E modulus of elasticity, lb/sq ft 

[F] structural influence coefficient matrix, streamwise angular 

deflection in radians at one station due to an applied load 
in pounds at another station 

(An element in the ith row and jth column indicates the 
angle of attack at station i due to a load at station j . ) 

G dimensionless load coefficient from leading edge of panel to an 

arbitrary point, — 

c av u 

AG dimensionless load coefficient over an arbitrary chordwise 

interval (eq. (9)) 

[I] spanwise integrating matrix for numerically integrating load 

coefficients 

M Mach number 

S total wing area including portion blanketed by fuselage, sq ft 

Sp area of one exposed wing panel, sq ft 

U free-stream velocity, ft/sec 

b wing span, ft 

c av average chord based on total wing area, ft 

c.p. chordwise center of pressure in local wing chords 

k constant obtained from series used for calculating rolling- 

moment or lift coefficients for a flexible surface (eq. (20)) 

m cotangent of sweepback angle of leading edge of panel 

p rolling velocity, radians/sec 
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A 

E 

F 
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dynamic pressure, lb/sq ft 
maximum thickness, ft 

longitudinal and lateral coordinates (fig. l) 
coordinates of apex of angle-of-attack panel (fig. l) 
coordinates of intersection of leading edge of panel and wing tip 
angle of attack, radians 

n/M 2 -1 

circulation, sq ft/sec 

total aileron deflection, measured in a plane perpendicular to 
the y axis, radians 

lateral coordinate in wing semi spans 

sweep angle, deg 

longitudinal coordinate in average: wing chords 

p(y - yj 

X “ X-l 

P(y - yt) 

x - x t 

rectangular matrix 
column matrix 
row matrix 

row matrix in which all elements ire unity 

Subscripts 


antisymmetric 

incremental value due to flexibility 
flexible case equilibrium value 
rigid case value 
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s 

symmetric 

le 

leading edge 

te 

trailing edge 

tip 

tip aileron 


ANALYSIS 


This section is concerned with the method of determining the result- 
ant aerodynamic force on a flexible wing with supersonic leading and 
trailing edges and streamwise tips due to such initial angle -of -at tack 
distributions as deflected ailerons and rolling velocities. First, the 
manner of obtaining a set of aerodynamic influence coefficients (defined 
as the loading coefficient at one station due to a unit angle of attack 
at another station) is given. From these coefficients, the aerodynamic 
loading and the resultant aerodynamic force due to an arbitrary distribu- 
tion of angle of attack in both the spanwise and chordwise directions may 
be obtained. Then, a set of structural influence coefficients, which 
will be assumed already available in the form of an angle of attack at 
one station due to a load at another station, will be combined with the 
aerodynamic influence coefficients to determine an equilibrium angle -of - 
attack distribution and resulting aerodynamic force. 


Aerodynamic Influence Coefficients 


The calculation of the aerodynamic influence coefficients is based 
on linearized potential flow theory. First, the equations are determined 
for the aerodynamic load coefficients at various regions of the wing due 
to a part of the wing, or panel, deflected at 1-radian angle of attack 
(fig. l). The panel has a supersonic leading edge which extends outboard 
at slope m from an apex at an arbitrary point x 1 ,y 1 to the wing tip. 
Networks of control points which correspond to load coefficient locations 
and to locations of apexes of the deflected panels are selected and a set 
of load coefficients for each panel location is calculated. These load 
coefficients are then converted to aerodynamic influence coefficient form 
by suitable superposition of the panels. This superposition is determined 
by the coefficients required to numerically differentiate an arbitrary 
angle -of -attack distribution. 

The equations used for obtaining the load coefficients in the various 
regions indicated in figure 1 due to a panel at a unit angle of attack are 
summarized here, while a more detailed presentation is given in appendix A 
The equations are first presented in the form of load coefficients which 
are proportional to the streamwise integration of pressures from the 
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leading edge of the panel to the desired lofd point, x,y. While somewhat 
simpler equations could have been written ir, terms of pressure coefficients, 
a brief numerical analysis indicated that a considerably larger number of 
stations would then have to be used to obtain a given accuracy. This was 
particularly true for cases with even moderately swept leading edges. 

The load coefficients for points x,y, which are located in regions 
affected only by the inboard portion of a panel at a unit angle of attack 
whose apex is at x y y x , can be expressed by the following equations: 

Region I (t < -l) 

G = 0 (!) 


Region II (-1 < t < l) 


G = 


2(x 




Pm 


rtp c av _s/p 2 m 2 -1 


-I A - 

cos 


+ T cosh 


-1 


1 

T 


( 2 ) 


where 0 < cos 1 ( ) < it. 
Region III (t > l) 


2(x - x 1 )(|3m - t) 
G — — * 11 ' * ■ — 

P c aWP 2m2 - ^ 


(3) 


These equations are valid for the three regions isolated from the wing 
tip (t+ < -l) and may be used to calculate loading coefficients for 
arbitrary angle -of -at tack distributions for plan forms whose edges are 
a.1 1 supersonic . 

For the remaining regions indicated by figure 1, the loss in loading 
due to the side edge must be taken into account. The following equation, 
expressed in terms of the distance from the tip of the leading edge of the 
panel x-^,y-t;, represents the loss in loading coefficients for 

Regions IV and V ("i 1 T t A) 


G = 


2(x - x-, ) f , -i 

== = A ( |3m - T t )cos 


rtP c aWP 2m2 " 1 


pm+ Tt(2pm + l) 

T t " P m 


+ 


a/pm(-T t )(l + T t )(pm+1) 


A) 


For region VI, which is defined by the equation 


x-x x > p(2y t - y x - y) 


(5) 



a separate load coefficient equation is used. For this region, it is 
convenient to combine the equation for the loading due to the inboard 
portion of the panel with the equation for the tip effect to give 
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G 


2 (yt -yJ 


:tc 


av 


7 y_y 7 

, -1 

2 (yt _ yi) , 

+ 2 

N 

L y-Vx 

\yt - yj 

bosh 

y-y x - 1 

I 1 yt - yi_ 


(6) 


The equations that give the load coefficients for a point x,y, which may- 
be located in various regions with respect to a panel at a unit angle of 
attack with an apex at x^y-^, are summarized in figure 1 . 

The load coefficient due to a panel at a unit angle of attack on the 
right-hand side of the wing center line can be calculated from the previ- 
ous equations for all points on both sides of the wing center line which 
are not affected by the wing tip on the opposite side from the deflected 
panel. For points which are affected by the opposite wing tip, an addi- 
tional correction would be required. However, for cases where this 
additional tip effect occurs only for points near the tip of the opposite 
side of the wing from the deflected panel, the loss in loading will be 
small and can often be neglected, and hence, equation ( 2 ) can still be 
used to obtain the load coefficient. The equation for this region is 
included in appendix A in order to determine limiting cases, however. 

From the calculated load coefficients for all load stations on the 
wing which result from a set of panel stations on one side, the m atrix 
for the antisymmetric load coefficients [G/J and the matrix of symmetric 
load coefficients [Gs] can be formed as follows: 

[Ga] = [G(x,y,x 1 ,y 1 )] - [G(x,-y,x 1 ,y 1 )] ( 7 ) 


[Gs] = [GCx^x^yj,)] + [GCxj-y^-^y-J] (8) 

These load coefficients represent loads resulting from the integration of 
pressures from the leading edge of a panel to a given load point, x,y. 

In order to obtain a chordwise distribution of load coefficient that 
can be used in conjunction with the structural influence coefficients, 
the operation 




is performed for each chordwise interval by subtracting the value of 
load coefficient at the back of each interval from the value at the 
front of the interval,. If desired, the differentiation could be calcu- 
lated more accurately through use of polynomial curve fitting (e.g. , 
ref. 12). Note that while the reduction of the G load coefficients to 
a chordwise incremental form, AG, corresponds to expressing them as 
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average pressures over the interval, a smoothing effect has been intro- 
duced. to make the subsequent numerical integration more accurate than if 
pressure coefficients had been used. The operation may be expressed 
symbolically by a matrix [B] as 

[AG] = [B][G] (9) 

Each column of the matrix [AG] represents tie distributed loading result 
ing from a panel such as an aileron at a unit angle of attack. The load 
coefficients for control surfaces vhich do rot coincide with the panel 
control points can be determined by interpolation of the panel values 
since the variations in these values are generally quite smooth. 

The AG panel load coefficients will now be adapted to the deter- 
mination of the loading due to a continuous arbitrary angle -of -attack ^ 
distribution. The following equation, in which the dimensionless vari- 
ables £,t) are used in place of x,y, will be evaluated numerically for 
the case of zero angle of attack at the plane of symmetry. 


■vte 


ag*U,t) = 




J o u le 






.(T)x) + 


- 


d -Si dT lx ( 1Q ) 


where AG*U,q) is the load coefficient distribution resulting from an 
arbitrary angle- of- attack distribution, a(fi,Tj x ). In order for the 
AGU^Ix^i) fu nction in integral form to correspond to the previously 
given AG in matrix form, a chordwise inte rval in | would have to be 
specified. This need not be done explicitly, however, since the equation 
is used only to illustrate the operations ]>erformed by the matrices. 

In evaluating equation (10) numerically, it will be convenient to 
define a quantity representing the product of the following terms and 
operations. 


A(E>TbSx>nx) 


AG(E,q,| 1 ,T] 1 ) 




+ 


K1 

asj 




1 


This quantity A, when expressed in terms of discrete intervals of the 
coordinates, is in the form of a load coefficient at one station due to 
a unit angle of attack at another station and will be called the aero- 
dynamic influence coefficient matrix. 

[A] = [AG][D] (H) 


The load coefficient due to an arbitrary angle -of -attack distribution can 
then be evaluated by use of a single matrix multiplication. 

{ AG*f = [A]{a} 
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The matrix [ D ] performs the operation t- — 


L 


( ) 


le 


a( ) 
as 


1J 


ty 


numerical differentiation and multiplies by the corresponding integration 
intervals d^dr^. The numerical differentiation coefficients used in 
this matrix are based on coefficients obtained by polynomial curve fitting 
(e.g., ref. 12) of the angle -of -attack curves. The degree of the polyno- 
mial employed depends upon the shape of the angle -of -attack curves and 
the accuracy desired. Each row of the matrix [D] represents the product 
of the integration interval and the coefficients which determine the rate 
of change of angle of attack that results from an arbitrary angle-of- 
attack distribution. Each column represents the product of the integra- 
tion interval and the coefficients of the slopes at the neighboring 
stations which produce a unit angle of attack at one station and zero 
angle of attack at the remaining stations. Thus, the aerodynamic influence 
coefficient matrix [A] represents the loading which results from a super- 
position of panels that give a unit angle of attack at one station and 
zero angle of attack at the remaining stations. 


A program for calculating the aerodynamic influence coefficients was 
set up for an IBM 650 digital computing machine. 


Structural Influence Coefficients 


For the analysis of this report, the structural influence coefficient 
matrix [F] will be assumed available from either experiment or theory in 
the form of a change in angle of attack at one station due to a concen- 
trated load at another station. Where possible, of course, it is desir- 
able to measure the coefficients directly from the actual structure. 

While a number of prediction methods are available (e.g., refs. 6 through 
10), the accuracy depends upon the complexity of the structure and the 
plan form. Although prediction methods which will yield angle -of -attack 
distributions due to applied loads are adequate for higher aspect ratios, 
they may not be adequate for low-aspect-ratio wings where chordwise bend- 
ing and root conditions become of greater importance. 


Calculation of the Stability Derivatives for a Flexible Wing 


In order to determine a particular stability derivative, the sets 
of aerodynamic and structural influence coefficients which were obtained 
for corresponding points on a plan form must be combined, and the equi- 
librium angles of attack and the resulting forces and moments on the 
flexible wing must be calculated. The equation for the incremental angle 
of attack due to flexibility, ccg, which results from an initial arbitrary 
distribution of angle of attack for the undeflected or rigid wing, op, 
may be expressed in integral equation form as 
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a£(t 1 ,n 1 ) - I s jf^ {wsi.nx.s.n)]^ 1 ^ 


te ao(e,n,ii»ni) 

Si 


*lx [ 3f ’ le 


<*»> 


SagUi>ni) 


■jd^d^j di dr, (12) 


where ap = ap + ap. 

When the previously discussed aerodynamic and structural influence 
coefficients are used, and a spanwise integrating matrix [i] is intro 
duced for numerically integrating the load coefficients in the final 
integration, equation (12) can be expressed in matrix form as 

)a E } = qS[F][l][A A ]{a K + a E f 0-3) 

The relaxation process described in reference 11, in which the incremental 
angle of attack |a E } is expressed in series form, is used to evaluate this 
equation. Once the incremental angle-of-attack distribution is obtained, 
the resulting lift and moment coefficients due to the rigid and incremental 
angles of attack can then be expressed in ai analogous series form. The 
antisymmetric case will be considered first. The rolling-moment coeffi- 
cient due to a rigid angle-of-attack distribution is given by the equation 

C Zr = LnJ[l][A A ])*R} ( l4 ) 

Since the incremental angle-of-attack distribution due to flexibility is 
small in comparison with the rigid value, a first approximation to the 
incremental angle of attack {a*;} is obtained through use of equation (13) 


= s[f][i][a a ]|»b( ( 15 ) 


For convenience, the incremental rolling-mcment coefficient is defined in 
terms of a unit value of dynamic pressure. 

C U = hJ[l][A A ]^} (16) 

The loading represented by causes an additional increment in the 

angle of attack due to flexibility which is determined by 

{P}- S [f][i)Ua{^} (17) 


and. 


= hJ[l][A A ]{pj 


(18) 
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The procedure is repeated so that the rolling-moment coefficient for the 
flexible case can then be expressed in series form as 

C l F = c Z R +Cz i q+C Z2 q. 2 + C Z3 q3+ . . . (19) 

After the relaxation process is repeated a sufficient number of times, 
the ratio of succeeding coefficients, Cj n /C2 n -i> approaches a constant, 
-k. The effect of succeeding terms can then be expressed in analytical 
form as C^ n q/(l+kq). For instance, if the value of Cj /c z and the 
succeeding ratios is a constant, the equation for rolling-moment coeffi- 
cient can be expressed as 


C Z F 


[, R) 1 + c Zr q 


+ 


C lz 

c h 


q 1 + 


C i3 <l/Cz 

1 + kq 



( 20 ) 


An analogous expression for lift coefficient can be obtained in a similar 
manner through use of the symmetric aerodynamic influence coefficients 
and the following equation. 


C L = LlJ[l][As]{af (21) 

The value of q = -l/k is the dominant characteristic or eigenvalue 
of dynamic pressure of the system of simultaneous equations (eq. 13). 

For an unswept or sweptforward wing, l/k is usually negative and repre- 
sents the dynamic pressure for which divergence occurs. For a sweptback 
wing, l/k is usually positive, and hence no real divergence speed is 
encountered. 


NUMERICAL EXAMPLE 
Description of Example 


Geometry . - The preceding analysis has been applied to the hypothetical 
supersonic configuration shown in figure 2. Pertinent geometric charac- 
teristics are shown in the figure. Since wing flexibility may have a 
major effect on lateral control, the following controls were investigated: 
a 0.143-semispan full-chord tip aileron, a 0.429-semispan l/3-chord 
trailing-edge aileron, and a slot -deflector spoiler with a chordwise 
center of pressure assumed to be acting along the 70-percent chord line 
and extending from 0.500 semispan to 0.857 semispan. 

Aerodynamic influence coefficients . - Aerodynamic influence coeffi- 
cients were calculated for two Mach numbers, 1.25 and 2.60, which are 
about the limits for which the lifting- surface theory should be useful 
for this plan form. The lower Mach number, for which the Mach wave sweep 
angle is 36.8° and the leading edge is only slightly supersonic, represents 
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a condition of large tip effect and interac :ions among the deflected q 
surfaces. The higher Mach number, with a Mich wave sweep angle of 67.4 , 
represents a condition of only small interaction effects. 

While a more complete analysis would have included an approximation 
of the effects of wing-body interference, these effects were neglected in 
this case since the principal emphasis was on lateral derivatives for 
which wing-body interference is small. The wing-body combination was 
replaced by a wing alone extending over the wing area blanketed by the 
fuselage. The only portion considered flexible was the part outboard of 
the wing-body juncture. However, the aerodynamic interaction effects 
between the flexible and rigid portions of the surfaces were still 
included. 

The exposed plan form was divided into a grid of 15 sections with 
5 equidistant spanwise intervals and 3 equidistant chordwise intervals 
(fig. 2 ). Since the aerodynamic loading varies considerably near the 
tip but is zero at the tip, the farthest outboard load station was a 
half interval from the tip. In addition, Iwo spanwise intervals were 
used for the portion of the wing blanketed by the fuselage. After the 
coordinates for the grid were selected, the load coefficients [GaI an< l 
[Gs] were calculated from equations (l) through (8) for both load and 
deflection locations. The matrix [B] for converting the [G] coefficients 
to the increment in load coefficient acting over each interval between 
adjacent chordwise grid points , and the mavrix [D] for differentiating 
the angle-of -attack distribution are given in tables II and III, respec- 
tively. In determining the matrix [D] , the numerical differentiation 
was calculated in the streamwise and constant percent chord directions 
since the wing stations for the influence coefficients were arranged in 
these directions. An examination of the structural influence coei^icients 
showed that the variations in angle of attack in the chordwise direction 
were less than those in the spanwise direction. Hence, a linear curve 
fit of two neighboring points was used to obtain coefficients for angle- 
of -attack slopes in the chordwise direction, while in the spanwise direc- 
tion, coefficients for slopes based on a pirabolic curve fit of three 
neighboring points (e.g. , ref. 12) were us 3d. Since only two neighboring 
points were used to determine the slopes ii the chordwise direction, the 
midpoints of the intervals were selected a 3 reference locations for the 
slopes. A more detailed description is gi/en in appendix B. The aero- 
dynamic influence coefficients were then calculated from equations (9) 
and (ll). 

Structural influence coefficients .- The structural influence coef- 
ficients, table I, were estimated from available measured influence 
coefficients of thin structures with somewhat similar plan forms. They 
are listed in a dimensionless form that is based on the deflection equa- 
tion from plate theory which is applicable to the relatively thin wings 
considered. These coefficients represent typical variations in deflec- 
tions of a homogeneous structure wing with the maximum thickness near the 
midchord. The addition of lateral controls to the structure was assumed 
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to be such that the structural influence coefficients were unchanged. 

The coordinates for the 15 deflection stations were located at 5 equally 
spaced intervals from the wing-body juncture to the wing tip and at the 
midpoints of 3 equal intervals in the chordwise direction (fig. 2). 

Calculation of the stability derivatives .- The rolling effectiveness 
of the trailing-edge aileron, tip aileron, and spoiler, and the damping 
in roll (C* p ), as well as the lift effectiveness (Ci^), were evaluated 
for a range of dynamic pressures for both Mach numbers through use of 
equations (13) through (21). The determination of the spanwise integrat- 
ing matrix [i] which appears in these equations is described in appendix B 
and listed in table IV. Calculations were also made using simplifications 
of the aerodynamic and structural characteristics described in the next 
section. 

In general, the convergence of the iteration process described in 
the analysis section was good in the calculation of results from this 
example. For the high Mach number, only two iterations of the rolling- 
moment or lift coefficient ratios described in the analysis section were 
needed to define the curves over the range of dynamic pressures shown. 

This was the same number required for the high-aspect-ratio configuration 
described in reference 1. For the low Mach number, the number of itera- 
tions required depended upon the rigid-wing load distribution. For the 
trailing-edge aileron, only two iterations were required, while for the 
other quantities , additional iterations were needed up to a maximum of 
five for C? when the complete aerodynamic and structural influence 
coefficients^were used. For all cases, the eigenvalue for the symmetric 
loading was essentially the same as that for the antisymmetric loading. 


Results and Discussion 


The calculations made of the lateral-control effectiveness and the 
stability derivatives Ci and Cp & over a range of dimensionless dynamic 
pressures for Mach numbers of 1.25 and 2.60are discussed and comparisons 
are made of the results by the complete method with those obtained by 
several aerodynamic and structural simplifications. 

Complete method results .- The ratios of stability derivatives for 
the flexible and rigid wings are shown in figure 3- The rigid values are 
listed in the following table to provide an additional comparison: 

M C La ^p ^te ^tip 

1.25 4.51 0.405 0.0981 o.o4n 

2.60 1.62 .197 .0334 .0247 

Tiie rigid values shown in the table were obtained by the numerical inte- 
gration and were in good agreement with corresponding analytical values. 
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This comparison indicates that the number of spamd.se stations selected 
vas adequate for this example. Since the chordwise locations of centers 
of pressure have considerable influence on the aeroelastic effects, the 
centers of pressure of the rigid-wing loadings are also shown (fig. )• 

The centers of pressure near the tip are of greatest interest since the 
largest structural deflections occur in this region* 

Looking first at the curves in figure 3 for the Mach number of 1 * 25 , 
the variations of flexible-rigid ratios with dynamic pressure are con- 
siderably different for the various derivatives. The trailing-edge 
aileron, with a rigid-wing center of pressure at the 0.83 chord, has the 
greatest loss with dynamic pressure. The loss for the spoilers with. a 
center of pressure at 0.70 chord, is considerably less, and the dynamic 
pressure for which control reversal occurs is approximately double tha 
of the trailing-edge aileron. The quantities whose centers of pressure 
for the rigid case are farthest forward (fig. *) , tip aileron Ci b , 

and C 7 decrease only slightly or even increase in effectiveness. An 
examination of the structural influence coefficients shows that the 
average rotation about an axis perpendicular to the plane of symmetry at 
each spanwise station may be represented api roximately by an elastic axis 
from about 0.30 chord near the tip to about 0.50 chord toward the root. 

The changes in rolling-moment coefficient die to flexibility are greatest 
when the differences between the approximate elastic axis location an 
the centers of pressure of the rigid loadings are the largest. 

The increase in with dynamic pressure shown in figure 3 (&) 

could be quite detrimental to the rolling performance of configurations 
with rolling-control effectiveness that decreases with dynamic pressure. 
The factors that contribute to this variation in C^ p can be explained, 
if they are examined in a manner similar to the series development use 
in the calculations. Since the center of pressure of. the rigid C Zp 
loading is fairly close to the approximate elastic axis, the average 
incremental angle of attack due to this loading is quite small. However, 
a significant amount of chordwise bending ii present and the center of 
pressure of the incremental load due to this bending is near the .7 . 

chord because of the relatively flexible afuerportion of the wing. This 
increment in angle of attack represents practically the total deflection 
mode shape at very low dynamic pressures, aid the resulting increment of 
loading will tend to decrease the magnitude of the rolling moment. The. 
mode shape at higher dynamic pressures can. oe visualized if the total wing 
distortion is separated into (l) this initial mode shape,. for w ic o y 
the magnitude will vary with dynamic pressure and which will produce a 
rolling-moment increment of opposite sign from the rigid value, and (2) 
the additional deflections due to loads and deflections that result from 
the initial mode shape and which produce a rolling-moment increment of 
the same sign as the rigid value. Since the center of pres sure. of the 
initial mode shape is approximately the sane as that of the spoiler, 
additional deflections will yield an incremental rolling moment which is 
of equal but opposite magnitude to that due to the initial mode shape a 
a dynamic pressure approximately the same as that for which rolling momen 
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reversal occurs for the spoiler. Therefore, the C rolling-moment 
ratio will increase to unity at this dynamic pressure. At higher values 
of dynamic pressure, the rolling -moment ratio for C^ p will be greater 
than unity. 

At the higher Mach number, changes with dynamic pressure would be 
expected to be similar but less because of the reduction in two-dimensional 
lift-curve slope. This trend is shown by the trailing-edge aileron and 
spoiler curves whose rigid-wing centers of pressure are the same for both 
Mach numbers. However, the values of Cp^, C p , and tip aileron Cjg, 
r ~\ 1 decrease with dynamic pressure because the centers of pressure of 
the rigid loadings are farther aft than for the low Mach number. 

Aerodynamic simplifications .- In order to determine under what 
conditions simplifications of the lifting-surface theory could be made , 
the following two simplifications were investigated. One simplification 
used was to assume that the magnitude of the load, at each station depended 
only on the local angle of attack and was equal to the loading obtained 
at that station if the entire wing was at the same angle of attack. Thus 
the aerodynamic influence coefficient matrix was a diagonal matrix with 
values equal to those for the wing at a unit angle of attack. The results 
obtained from these values will be referred to as modified strip theory 
results. An additional simplification of the modified theory was made 
by letting the three aerodynamic coefficients at each spanwise station be 
of equal value, but leaving their sum the same as in the previous case. 
Thus, the chordwise distribution of loading would be proportional to the 
local angle of attack, which is analogous to the two-dimensional case. 
Results obtained from this additional modification will be referred to as 
simplified strip theory results. 

Comparisons of the flexible-rigid ratios from the modified strip 
theories with those from the complete aerodynamic theory are shown in 
figure 5„ At the lower Mach number, the values of rolling-moment and 
lift coefficient ratios from the strip theories are seen to be much less 
than those from the lifting-surface theory. This discrepancy is par- 
ticularly true for the simplified strip theory case which underestimates 
the values of Cp p by approximately 50 percent at the higher dynamic 
pressures. As mentioned previously, the effect of the tip as predicted 
by lifting-surface theory is to reduce the loading near the tip, with 
the greatest reductions occurring over the afterportions. Hence, if a 
strip theory based on a constant angle-of -attack distribution is used, 
the results will underestimate tip effects due to wing distortions since 
these angle-of -attack distributions are concentrated more toward the tip. 
Thus the incremental rolling moments due to flexibility predicted by the 
strip theories will tend to be larger so that the resultant rolling- 
moment ratio will decrease more with dynamic pressure. Another reason 
for the deviations in is that the chordwise center of pressure of 

the rigid loading from the strip theories is farther aft than that from 
lifting- surface theory and therefore it is farther from the approximate 
elastic axis. 
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The modified strip theory results for show good agreement with 

the lifting-surface theory even for the lower Mach number, since the rigid' 
wing loadings are the same for the two theories and the variation of 
C La F / C L aR vith dynamic pressure is small for this rigid- wing center-of- 

pressure location as it is close to the approximate elastic axis. How- 
ever, the center of pressure of the rigid Cp a loading predicted by the 
simplified strip theory is a little farther aft than that predicted by 
lifting-surface theory; therefore, a slightly greater reduction of 
CL a p/CL(x-^ with dynamic pressure occurs. For the trailing-edge aileron, 

predictions of reversal dynamic pressure obtained by use of the modified 
strip theories are about 25 percent lower than those predicted when 
lifting-surface theory is used. 

At the higher Mach number , good agreement is obtained between the 
strip theory and lifting-surface theory results because the tip and other 
aerodynamic interaction effects extend over only small spanwise portions 
of the wing. 

Structural simplifications . - To get an indication of the importance 
of chordwise bending for this case, a simplification of the structural 
matrix was made by assuming that each spanwi ;e section rotated in the 
streamwise direction as a unit. Then, each _oad was assumed to produce 
an angle of attack at a spanwise station equal in magnitude to the 
average of the values of the deflections at c.he three chordwise stations. 

Comparisons of the flexible -rigid ratios with and without chordwise 
bending are shown in figure 6. At both Mach numbers, neglecting chord- 
wise bending makes an appreciable difference in the aeroelastic effects. 
The aileron reversal dynamic pressure is underestimated by about 15 per- 
cent while decreases in Cj p and Cx^ also occur in both cases. The 
reason for this trend can be explained as foLlows: An examination of 

the structural influence coefficients (table i) indicates that the 
structure is relatively flexible over the afterportion of the wing in 
that deviations from the average deflection tend to be greater in the 
afterregions and also loads applied to the afterportions produce greater 
deflections. At low dynamic pressures, the average value of incremental 
angle of attack due to deflections will be the same for the two cases 
with and without chordwise bending, but the angle of attack will be 
concentrated farther aft for the chordwise tending case. These angles 
of attack will produce approximately the sane increment of rolling-moment 
coefficient which will tend to counteract tie originally applied load. 
However, as dynamic pressure is increased, these incremental loads will 
produce an additional angle -of -attack change. For the case where chord- 
wise bending is neglected, the latter angle- of— attack increment will be 
smaller t han the complete structure case since the loads are smaller over 
the afterportion of the wing. Consequently, the loading resulting from 
the latter angle of attack, which is in the direction of the originally 
applied loading, is smaller than the complete structure case. Thus the 
rolling-moment coefficient ratio will decrease more at the high dynamic 
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pressures when chordwise bending is neglected. The lifting-surface 
effects of Mach number are relatively small in this comparison since , as 
may be deteimined from figure 6, the values of the flexible -rigid ratios 
for the cases with and without chordwise bending are approximately the 
same for equal values of the parameter q/p, This parameter takes into 
account the differences in two-dimensional lift-curve slopes for the two 
Mach numbers. 

It is of interest to note that, when chordwise bending is included, 
the magnitude of the constant k is greater since the corresponding 
chordwise angle -of -attack distribution and load distribution are farther 
aft and hence a larger deflection and load occurs. Thus the resulting 
contribution of the chordwise bending in this case is to reduce the 
possibility of a divergent condition. For a case with the largest deflec- 
tions toward the forward part of the wing, a trend toward divergence 
would be expected. 


CONCLUDING REMARKS 


A method has been presented for calculating lifting -surface aero- 
dynamic influence coefficients for swept wings with supersonic leading 
edges. Through use of these aerodynamic influence coefficients, together 
with a set of structural influence coefficients, aeroelastic effects were 
calculated for an example low-aspect-ratio wing for Mach numbers of 1.25 
and 2,6 o. An evaluation of results from these example calculations for 
several lateral controls and stability derivatives has led to the 
following conclusions: 

1, Comparisons of results from a modified strip theory with those 
from the lifting- surface theory indicate that the modified strip theory 
was satisfactory for the high Mach number case, but inadequate for the 
low Mach number case. A large part of the discrepancy was due to the fact 
that the modified strip theory did not predict with sufficient accuracy 
the chordwise distribution of loading near the tip which resulted from 
the wing deflections. 

2. For this example, neglecting chordwise bending tended to predict 
greater reductions in rolling-moment and lift coefficients with dynamic 
pressure than when the complete structural characteristics were used. 

3* The trailing-edge aileron had the greatest loss in effectiveness 
of all the controls studied and resulted in a relatively low reversal 
speed. The loss in lateral-control effectiveness of a slot-deflector 
spoiler was considerably less, having a reversal dynamic pressure 
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approximately double that of the trailing-eclge aileron. The tip aileron 
had the least variation in control effectiveness with dynamic pressure, 
and at the lower Mach number the effectiveness actually increased at 
high dynamic pressures. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Jan. 21, 1959 
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APPENDIX A 

DERIVATION OF EQUATIONS FOR AERODYNAMIC LOADING 


The equations for the loading due to the wing panel at 1-radian angle 
of attack shown in figure 1 will be developed. This development is 
equivalent to the determination of the loading due to an aileron deflec- 
tion for various regions over a wing including the regions of interaction 
effects between the inboard end of the aileron and the tip. Linearized 
potential flow theory is used. While , of course, part of this develop- 
ment has been given in the literature previously, the derivation of all 
the equations is indicated for completeness. In determining the result- 
ing loads due to an arbitrary angle -of' -attack distribution, it was found 
desirable to express the aerodynamic coefficients in terms of circulation 
coefficients rather than pressure coefficients, since the total load could 
then be obtained with greater accuracy provided a sufficient number of 
control points were used to define the angle-of -attack distribution. This 
was particularly true for cases with swept leading edges. 

The equations for the pressure on the top surface for the inboard 
portion of a panel at 1-radian angle of attack were given in reference 13 
and are expressed by the notation shown in figure 1 as 

pC p =0 (r < -1) (Al) 


PC p 


2pm -l 1 - PmT 

— : cos 

Wp 2 m 2 - 1 Pm - t 


(-1 < T < l) 


(A2) 


where 0 < cos 


1 - $mr \ 
Pm - t J 


< 


Jto 


PC T 


2pm 


•J p 2 m 2 - 1 


(t > 1) 


(A3) 


The last equation represents the pressure over the portion of the panel 
which is undisturbed by either edge. 


To get the dimensionless load coefficient G = r/Uc av between the 
leading edge of the panel and a point (x-xj.), the following equation 
was used to integrate the previous equations in the streamwise direction. 



iX - X! 


(x - Xi) 


C-nd(x- Xi) 


le 


(A4) 
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The equations for the load coefficient become 

G = 0 (t < -l) (A5) 


G = 


2(x ~ xi) 

It P c av 


gm - r 


= cos 

1 


-l 


1 - pmr 
Pm - r 


+ T cosh -1 



(-1 < T < 1) 

(A 6) 


2(x - x x ) (pm - t) 

P c av *iP 2 m 2 “ 1 


(t > l) 


(A7) 


Equations (A5), (A6), and (A7) can be used to calculate the load coeffi 
cient distribution due to an arbitrary angle ~of -attack distribution over 
plan forms whose edges are all supersonic. 

In order to obtain the loading over the portions of the panel near 
the tip, an additional term is needed which will be expressed as a correc- 
tion to the load coefficient previously obtained for the inboard section 
of the panel. The method of Eward (ref. l4) was used in which the . 
effect of the subsonic edge is obtained by a cancellation of a portion of 
the integration over the wing surface in the equation for the velocity 
potential. The dimensionless load coefficient G(x,y), which is propor- 
tional to the loading from the leading edge to the point x,y, is defined 
in terms of the velocity potential and circulation function as follows: 


G = 



2cp 


Uc 


av 


(A8) 


For regions where no interactions exist between the Mach line 
reflected from the tip and the inboard end cf the panels (regions IV and 
V in fig. l), the area of integration is shewn in the following sketch. 
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The loss in load coefficient due to the tip was expressed by evaluating 
the load coefficient obtained from integrating over the shaded area shown 
in the sketch and subtracting from it a value obtained for the panel 
without the free edge. 


G = 0 (r t < -1) 


(A9) 


2(x - xj 


rtPeav^P 2 ™ 2 ~ 1 


-(pm- T-j- ) cos 


-l 


Pm+ t.(- (2pm + 1) 

T t " P m 


2 s/pm(-T t )(l + T t )(pm+l)j- (-1 < T t < 0 ) 


(A10) 


where 0 < cos“ ( ) < at. 


An additional modification of the tip effect must be made for cases 
where the inboard edge of the panel is relatively close to the tip so 
that the reflected Mach line from the tip intersects the inboard edge. 
This condition is given by the equation 


X “Xi > p(2y t -y 1 -y) (All) 

The region of integration is shown in the following sketch. An examina- 
tion of the sketch shows that the values of load coefficient for a 



point x,y in region VI of figure 1 are the same as the load coefficient 
at the boundary between regions V and VI for the same value of y. The 
load coefficient from the leading edge to a point in region VI can be 
obtained by substituting the value of the boundary 


X-Xi = p(2y t -y x -y) 
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into equations (a 6) and (A10). These equations can then he combined to 
give 


y- y x 
7t - 


-i 2 (yt~yi 

y- y x 


- 1 + 21 - 


yt " yi 


Note that this equation is independent of the leading-edge sweep angle 
as may be seen from the area of integration shown in the previous sketch. 

The equations for the various regions shown in figure 1 are summarized 
in the main body of the report. 

For cases with very low aspect ratios and leading edges just slightly 
supersonic, another equation must be developed for points on the opposite 
side of the wing panel from the deflected panel which are affected by the 
panel and which are defined by the equation 

x -xx > p(2yt + y x + y) (A!3) 

The area of integration is indicated by the shaded region in the follow- 
ing sketch. 





The load coefficient obtained for this regi Dn is 
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G = 


JtC 


av 




m x + y 

jrr 


■nr 


cos 


-i 


~mx -y + 2pm(y-y t ) +2y t ' 


mx + y 


^(yt-y)I"*My^i-^t),,!nIl 1/2 . (y - Jl)cosh -i 

Pm + 1 J 


2(yt -yi) 


y t -y 


r-‘i) 


(A14) 


Fortunately, for most cases, including the lowest Mach number case con- 
sidered for the example used in this report, the effect of the opposite 
wing tip from the deflected panel is small, and the load coefficient for 
this region may be obtained by using equation (A6). This approximation 
results in including the effect of the extra triangular area near the 
apex of the panel shown in the previous sketch. 
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APPENDIX B 

DIFFERENTIATING AND INTEGRATE FACTORS USED 
FOR NUMERICAL EXAMPLE 


Two sets of numerical coefficients will be evaluated for the deter- 
mination of the loading over the flexible w .ng: (l) the differentiating 

matrix [D] for differentiating the angle of attack (eq. (ll)), and (2) 
the spanwise integrating matrix [i] for integrating the loading coeffi- 
cients (eqs. (13 ), (l4), and (2l)). 


Differentiating Matrix [D] 


This portion of the appendix will be concerned with the method used 
to convert the panel incremental load coef f Lcients to aerodynamic influence 
coefficients (eq. (ll)), that is, coefficients of the form of the incre- 
mental load coefficient at one station due :o a unit angle of attack at 
another station, the angle of attack of all other stations being zero. 

This will make the aerodynamic coefficients compatible with the struc- 
tural influence coefficients which are expressed as an angle -of -attack 
distribution. The angle of attack at one station will be defined in 
terms of the slopes of the neighboring points in both the chordwise and 
spanwise directions from which the resulting loading may be obtained by 
use of panel angle-of -attack coefficients. Polynomial curve-fitting 
methods were used to obtain the slopes (ref. 12), with the degree of the 
polynomial depending on the accuracy desire i. 

For the example case, the angle -of -attick variations in the chord- 
wise direction were much less than those in the spanwise direction so 
that simpler coefficients were used to express the variations in the 
chordwise direction. The chordwise variatiDns will be considered first. 
Three equally spaced intervals were used for the exposed portion of the 
wing panels (fig. 2). Since the [AG] coefficients represent loadings due 
to a panel at a unit angle of attack extending rearward, the arbitrary 
angle-of -attack distribution is first converted to an incremental angle- 
of-attack form starting from the leading ed^e (i.e., coefficients consist- 
ing of the slopes at each interval times ths interval are formed). For 
this case the increments can be expressed simply as 

at - cc. 


Denoting the leading edge of front, mid, an I rear intervals by the sub- 
scripts f, m, and r, respectively, the incremental angles of attack may 
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be written as 


Ax^ = a_£> 

+ a m > 
Ax r - -am + a r 


(Bl) 


The left-hand side of equations (Bl) denotes a chordwise row location in 
table III, while the right-hand side indicates a chordwise column loca- 
tion and gives the multiplying factor to be combined with the appropriate 
spanwise coefficient which will now be evaluated. 


For the spanwise variation, the eight control stations used for the 
angle -of -attack and slope locations are given in the following table. 


Station 

0 

1 

2 

3 

4 

5 

6 

7 

a at i) 

0 

o.i43 

0.286 

0.429 

0.572 

0.715 

0.858 

1.000 

da/ dT] at tj 

0 

• 143 

.286 

.429 

• 572 

• 715 

00 

Lf\ 

OO 

.929 


The outboard station for the slope was chosen at a half interval from the 
tip since the value at the tip would have no effect on the aerodynamic 
load. The first three angle-of -attack stations denote the portion of the 
wing blanketed by the fuselage and are used only for angle-of -attack 
distributions for the rigid case. Hence, the values given in table III 
apply only for cases where no chordwise variation in angle of attack 
occurs for these three stations. The slope at each station was obtained 
by passing a parabola through three adjacent spanwise stations and dif- 
ferentiating the result. The resulting coefficients consisted of the 
slope times the appropriate interval [ Atj (da/dT] ) ] , and are listed as 
follows : 

‘"N 

( At) = ! ! a ° + 2a i if) 





(B2) 
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The left-hand side of equations (B2) indicates a spanwise row location in 
table III, while the right-hand side denotes a spanwise column location 
and the corresponding spanwise coefficient ho be multiplied with a chord- 
wise coefficient given by equations (Bl). 


Integrating Matrix [i] 


The matrix [i] represents a spanwise integrating matrix for summing 
the load coefficients (eqs. (l3)j (l^-), and (2l)). Simpson* s rule factors 
(e.g., ref. 12) which are based on passing a parabola through three 
adjacent points, were used to integrate over the interval from rj = 0 
to rj = O.858. In order to better approximate the shape of the loading 
over the interval from r) = O.858 to the tip, an equation of the following 
form was used to obtain the integrating factors. 

G(T]) = k 0 (l- 7 ] ) 1/2 + k 1 (l-T ] ) 3/2 

The constants k 0 and k 1 were evaluated in terms of the value of G(t)) 
at T] = O.858 and 0.929* The integration constants were then determined 
from the above equation. The integration constants for the entire span, 
including the appropriate interval, are given in table IV. 
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TABLE I.- STRUCTURAL INFLUENCE COEFFICIENTS 
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Station location for structural 



Taper ratio = 0.4 aileron 

Body span to wing 
span ratio = 0.286 


Figure 2.- Geometric characteristics of numerical example. 
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(b) M = 2. So 

Figure 4 .- Centers of pressure of rigid-wing loadings 
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(b) M = 2.60 

Figure 5.- Comparisons of lifting surface theory results with strip 

theory results. 







Flexi ble Flexible 

Rigid _ Rigid 


'igure 6 



(a) M = 1,2 5 



o .01 .02 .03 .04 .05 

q/E (t/c ) 3 


(b) M = 2.60 

. - Effect of chordwise bending on rclling and lift effectivenes 


NASA - Langley Field, Va. A -1^9 




